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Abstract 

The problem addressed in this paper is to control a plant so as to have its output 
tracking (a family of) reference commands generated at a remote location and trans- 
mitted through a communication channel of finite capacity. The uncertainty due to 
the presence of the communication channel is counteracted by a suitable choice of the 
parameters of the regulator. 
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1 Introduction 

In distributed control systems, sensors, actuators and control unit may be placed at loca- 
tions which are geographically separated. Information among these devices must then be 
exchanged through a finite bandwidth channel. 

The problem addressed in this paper is to control a plant so as to have its output tracking 
(a family of) reference commands generated at a remote location and transmitted through 
a communication channel of finite capacity. What renders the problem in question different 
from a conventional tracking problem is that the tracking error, that is the difference be- 
tween the command input and the controlled output, is not available as a physical entity, as 
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it is defined as difference between two quantities residing at different (and possibly distant) 
physical locations. Therefore the tracking error as such cannot be used to drive a feedback 
controller, as it is the case in a standard tracking problem. 

The actual tracking error not being available, it is natural to approach the problem by 
reconstructing the tracking error starting from the information transmitted through the 
communication channel. For the reconstruction to successfully take place, the information 
must be suitably encoded. A possibility is to make use of a sufficiently large number of bits, 
so as to render the magnitude of the difference between the true reference signal and the 
reconstructed one negligible. However, in the framework of distributed control systems, the 
constraint on the available bandwidth is usually tight and adopting encoding schemes which 
require a large number of bits may not be practically feasible. The approach pursued in 
this paper is rather to counteract the uncertainty due to the presence of the communication 
channel by a suitable choice of the parameters of the regulator. 

In Section 2, the formulation of the problem is made more precise, whereas the procedure for 
encoding the reference command is described in Section 3. Section 4 introduces the regulator 
which guarantees the achievement of the control goal using the reconstructed tracking error. 
The main results of the paper are stated and proved in Section 5. The proof consists of two 
steps. First, boundedness of the closed-loop trajectories are shown (Section 5.1), and then 
asymptotic convergence to zero of the tracking error is concluded (Section 5.2). The results 
are illustrated by an example in Section 6. 



2 Problem statement 

Generally speaking, the problem in question can be defined in the following terms. Consider 
a single-input single-output nonlinear system modeled by equations of the form 

i = f{x)+g{x)u 
y = h{x) 

and suppose its output y is required to asymptotically track the output y^cs of a remotely 
located exosystem 

w = s{w) w G W , . 

2/dos = yriw). 

The problem is to design a control law of the form 

u = 9{^,y,w^) ^ ^ 

in which represents a sampled and quantized version of the remote exogenous input w, so 
as to have the tracking error 

eit)=yit)-y,iwit)) (4) 
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asymptotically converging to zero as time tends to oo. Note that the controller in question 
does not have access to e, which is not physically available, but only to the controlled output 
and to a sampled and quantized version of the remotely generated command. 

We will show in what follows how the theory of output tracking can be enhanced so as to 
address this interesting design problem. In particular, we will show how, by incorporating 
in the controller two (appropriate) internal models of the exogenous signals, the desired 
control goal can be achieved. One internal model is meant to asymptotically reproduce, at 
the location of the controlled plant, the behavior of the remote command input. The other 
internal model, as in any tracking scheme, is meant to generate the "feed-forward" input 
which keeps the tracking error identically at zero. 

We begin by describing, in the following section, the role of the first internal model. 

3 The encoder-decoder pair 

In order to overcome the limitation due to the finite capacity of the communication channel, 
the control structure proposed here has a decentralized structure consisting of two separate 
units: one unit, co- located with the command generator, consists of an enco(ier which extracts 
from the reference signal the data which are transmitted through the communication channel; 
the other unit, co-located with the controlled plant, consists of a (decoder which processes the 
encoded received information and of a regulator which generates appropriate control input. 

The problem at issue will be solved under a number of assumptions most of which are 
inherited by the literature of output regulation and/or control under quantization. The 
first assumption, which is a customary condition in the problem of output regulation, is 
formulated as follows. 

(AO) The vector field s(-) in Q is locally Lipschitz and the initial conditions for ^ are 
taken in a fixed compact invariant set Wq. < 

The next assumption is, on the contrary, newer and motivated by the specific problem 
addressed in this paper. In order to formulate rigorously the assumption in question, we 
need to introduce some notation. In particular let \x\s denote the distance at a point x G 
from a compact subset S C M", i.e. the number 

\x\q := max \x — y\ 

and let 

Lo = max \xi - yi\ . (5) 

2e[l,...,r] 
(x,y)(^WQXWQ 

Furthermore, having denoted by A'^f, the number of bits characterizing the communication 
channel constraint, let be the largest positive integer such that 

N,>r\\og^N] (6) 

where [w], w G M, denotes the lowest integer such that [w] > v. 
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With this notation in mind, the second assumption can be precisely formulated as follows. 
(Al) There exists a compact set W D Wq which is invariant for w = s{w) and such that 



\w\wo>Vr 



2N 



W being compact and s(-) being locally Lipschitz, it is readily seen that there exists a 
non decreasing and bounded function M(-) : M>o —>■ M>o, with M(0) = 1, such that for all 
wio e W and W2q G W and for alH > 

\wi{t) - W2{t)\ < M{t)\wio - W2o\ (7) 

where wi(t) and W2{t) denote the solutions of © at time t passing through wio and, respec- 
tively, 1^20 at time t = 0. 

This function, the sampling interval T, the number Lq defined in © and the number 
fulfilling 0, determine the parameters of the encoder-decoder pair, which are defined as 
follows (see (TU], |H], [SI for more details). 

Encoder dynamics. The encoder dynamics consist of a copy of the exosystem dynamics, 
whose state is updated at each sampling time kT, e N, and determines (depending on the 
actual state of the exosystem) the centroid of the quantization region, and of an additional 
discrete-time dynamics which determines the size of the quantization region. Specifically, 
the encoder is characterized by 

L{k) 



L{k+1] 



N ^ ^ 



We{kT) = We{kT ) -|- Wq(/c)- 

m = Lo 



N 



in which Wq represents the encoded information given by, for i = 1, . . . 



Wq,i{k) = sgn{wi{kT) - We,i{kT )) 



N\w,{kT) - w,,,{kT-)\ 

m 

N\wi{kT) - w,4kT-)\ 

m 



N even 



iVodd. 



At each sampling time kT, the vector Wq{k) is transmitted to the controlled plant through the 
communication channel and then used to update the state of the decoder unit as described in 
the following. To this regard note that each component of the vector Wq{k) can be described 
by [log2 N] bits and thus the communication channel constraint is fulfilled. 

Decoder dynamics The decoder dynamics is a replica of the encoder dynamics and it is given 
by 

Lik) 



L{k 



-I) 



r^^L[k) 



Wd{kT) = Wd{kT 

m = Lo 



Wq{k)'- 



N 



Wd{0-)=w,{0- 
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If, ideally, the communication channel does not introduce delays, it turns out that Wd{t) = 
Wc{t) for all t > 0. Furthermore, it can be proved that the set W characterized in Assumption 
(Al) is invariant for the encoder (decoder) dynamics and that the asymptotic behavior of 
We(t) {wd{t)) converges uniformly to the true exosystem state w(t), provided that T is 
properly chosen with respect to the number N and the function M(-). This is formalized in 
the next proposition (see [S], [S] for details). 

Proposition 1 Suppose Assumptions (AO)-(Al) hold and that the sampling time T and the 
number N satisfy 

N>y/¥M{T). (9) 

Then: 

(i) for any Wd(0~) G Wq and w(0) G Wq, Wd{t) G W for all t > 0; 

(a) for any Wd{0") G Wq and w{0) G Wq, 

lim \w{t) - Wd(t)| = 

with uniform convergence rate, namely for every e > there exists T* > such that for all 
initial states Wd(0~) G Wq, w{0) G Wq, and for all t > T* , \w(t) — Wd(t)\ < e. 

Proof. As W is an invariant set for w = s{w), the proof of the first item reduces to show 
that, for all k > 0, if w^ikT') G W then necessarily w^ikT) G W . For, note that this is 
true for /c = 0. As a matter of fact, since tyd(0~) G Wq C W and by bearing in mind the 
definition of w;q, it turns out that |u;d(0) — w{Q)\ < ^Lq/2N which implies, by definition of 
W in Assumption (Al), that ifd(O) G W . For a generic /c > note that, again by definition 
of Wq, it turns out that \wd{kT) — w{kT)\ < y/rL{k)/2N. But, by the second of (jH)) and by 
condition Q, L{k) < L{k — 1) < Lq yielding |wd(A;T) — w{kT)\ < ^JtLqIIN which implies 
Wd{kT) G W. This completes the proof of the first item. The second item has been proved 
in jH], 0. < 

Remark. By composing ® with Q it is easy to realize that the number of bits Nb and the 
sampling interval T are required to satisfy the constraint 

Nb>r\\og,{V^M{T))] (10) 

in order to have the encoder-decoder trajectories asymptotically converging to the exosystem 
trajectories. Since the function M(-) depends on the exosystem dynamics and on the set 
Wq of initial conditions for equation ()1U|) can be interpreted as a relation between the 
bit-rate of the communication channel and the exosystem dynamics which must be satisfied 
in order to remotely reconstruct the reference signal. 
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4 The regulator 



4.1 Basic hypotheses 

As in most of the literature on regulation of nonlinear system, we assume in what follows that 
the controlled plant has well defined relative degree and normal form. If this is the case and 
if the initial conditions of the plant are allowed to vary on a fixed (though arbitrarily large) 
compact set, there is no loss of generality in considering the case in which the controlled 
plant has relative degree 1 (see for instance P). We henceforth suppose that system ^ is 
expressed in the form 

z = f{z,y,ij) zeW 
y = qiz,y,fi) + u yeR 

in which /i is a vector of uncertain parameters ranging in a known compact set P. Initial 
conditions {z{0),y{0)) of are allowed to range on a fixed (but otherwise arbitrary) 
compact set Z X y C X M. 

It is well known that, if the regulation goal is achieved, in steady-state (i.e when the 
tracking error e{t) is identically zero) the controller must necessarily provide an input of the 
form 

Uss = LsVriw) - q{z,yr{w), fl) (12) 

(where Lsy^{-) stands for the derivative of ?/r(') along the vector field s(-)) in which w and z 
obey 

fi = 

w = s{w) (13) 

Z = /(2,l/r(w),/i) . 

As in |2], we assume in what follows that system (fT^ has a compact attractor, which 
is also locally exponentially stable. To express this assumption in a concise form, it is 
convenient to group the components fi,w,z of the state vector of (jl3j) into a single vector 
z = col(yU, w, z) and rewrite the latter as 

z = fo(z) . 

Consistently, the map (fT^ is rewritten as 

Uss = qo(z) , 

and it is set Z = P X ly X Z. The assumption in question is the following one ^ 

^Recall that, if the positive orbit of a compact set X of initial conditions of a system 

i = fix) (14) 

is bounded, the w-hmit set of X under the flow of (|14|l - denoted uj{X) - is a nonempty, compact invariant set 
which attracts X uniformly. If (-u{X) is in the interior of X, then X is also stable in the sense of Lyapunov. 
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(A2) there exists a compact subset Z of P x W x MP which contains the positive orbit of 
the set Z under the flow of (jl3p and uj{Z) is a differential submanifold (with boundary) of 
P X W X W^. Moreover there exists a number di > such that 

zePxWxW\ \z\^{z)<di zgZ. 

Finally, there exist m > 1, a > and d2 < di such that 

zoePxW xW, |zoU(z) <d2 =^ |zWL(z) < "^e-°*|zoU(z) , 

in which z(t) denotes the solution of (fT5|) passing through zq at time t = 0. < 

In what follows, the set w{Z) will be simply denoted as ^o- The final assumption is an 
assumption that allows us to construct an internal model of all inputs of the form Uss{t) = 
qo(z(t)), with z(t) solution of (fT^ with initial condition in ^o- This assumption, which can 
be referred to as assumption of immersion into a nonlinear uniformly observable system, is 
the following one. 

(A3) There exists an integer d > and a locally Lipschitz map : M'' ^ M such that, for 
all z E Aq, the solution z(t) of (jl3|) passing through zq at t = is such that the function 
u{t) = qo(z(t)) satisfies 

u^''\t)+ip{u^''-^\t),...,u^^\t),u{t)) = 0. < (15) 
4.2 The design of the regulator 

Using the technique described in j2] , the first step is the construction of an internal model for 
(fT^. viewed as an autonomous system with output (fT^ . To this end, consider the sequence 
of functions recursively defined as 

n(z) = qo(z), ri+i(z) = — ^fo(z) 

oz 

for i = 1, . . . ,d — 1, with d as introduced in assumption (A3), and consider the map 

r : P xW xW R'^ 

{H.w.z) col(ri(z),r2(z), . . . ,rd(z)) . 

If /c, the degree of continuous differentiability of the functions in flllj) . is large enough, the 
map r is well defined and C^. In particular t{Aq), the image of under r is a compact 
subset of W^, because Aq is a compact subset oi P xW xW^. 

Let (y^c : M"^ — M be any locally Lipschitz function of compact support which agrees 
on t{Aq) with the function defined in (A3), i.e. a function such that, for some compact 
superset S of t{Aq) satisfies 

Vc(^) = for all ?7 ^ 5 
Vci'n) = 'fiv) for all V e t{Ao). 
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With this in mind, consider the system 

e = $c(0 + G(Ms.-ro 



(16) 



in which 



$c(0 



6 
6 



V-V5c(6,6, 



G 



( K'Cd-l \ 



1 



the Cj's are such that the polynomial A'^ + cqA 



■ + Cd_i = is Hurwitz and k is a positive 



number. As shown in |2I, if k is large enough, the state ^(t) of ()16p asymptotically tracks 
T(z(t)), in which z(t) is the state of system (fTSj) . Therefore r^(t) asymptotically reproduces 
its output (|12j). i.e. the steady state control Uss{t). As a matter of fact, the following result 
holds. 



Lemma 1 Suppose assumptions (Al) and (A2) hold. Consider the triangular system 

z = fo(z) 

i = $,(0 + G(qo(z)-ro. 



(17) 



Let the initial conditions for z range in the set Z and let S be an arbitrarily large compact 
set of initial condition for ^. There is a number k* such that, if k > k* , the trajectories of 
are bounded and 



graph (r 



^0 



u Z X E) 



In particular gTa.ph.{T\j^^g) is a compact invariant set which uniformly attracts ZxS. Moreover 
graph(r|_4„) is also locally exponentially attractive. 

Note that, as a consequence of assumption (Al), of the fact that r(-) is a continuous 
vector field and of the fact that the compact set S can be taken arbitrarily large, it is 
possible to assume, without loss of generality, the existence of a positive ^2 such that 



z G P X X 



(z,OU(ZxH) < d2 



(z,OgZx 



In view of Lemma ^ it would be natural - if the true error variable e were available for 
feedback purposes - to choose for (jlip a control of the form 

u = Tt, — ke, 

with k a large number. This control, in fact, would solve the problem of output regulation 
(see 12]). The true error e not being available, we choose instead 

e = y-yriwd) (19) 
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and the controller accordingly as 



u = T^-ke. 

The result which will be proven next is that there exists k* > such that if /c > /c* the 
regulator designed above solves the problem in question (provided that and T satisfy the 
condition of Proposition [T)) . 



5 Main results 



5.1 Trajectories of the closed loop system are bounded 

To prove that the proposed regulator solves the problem, we show first of all the trajectories 
of the controlled system, namely those of the system 



Wd = s[Wd) 

y = qiz,y,iJ,) + - k{y -y,{wd)) 
e = MO-Gk{y-y,{wd)) 



Wd{kT) = Wd{kT-) + w^{k) 



m 

2N 



(21) 



are bounded. To study trajectories of (j2Tj) it is convenient to replace the coordinate y by 

e = y- VriWd) 

to obtain the system 



Wd = S{Wd) 

Z = /(2,e + ?/r(Wd),/i) 

e = MO-Gke 

e = q{z,e + y,{wd),fi) - Lsyriwd) + - ke. 



(22) 



This system can be further simplified by changing the state variable ^ into = ^ — Ge and 
setting p = col(/i, Wd, z,C,), so as to obtain a system of the form 



p = Fo{p) + Fi{p, e)e 

e = Ho{p) + Hi{p,e)e - ke , 



(23) 



in which 



Foip) 



\ 





s{wd) 

/(2,l/r(Wd),/i) 

V $c(0 + G{-q{z, y,{wd), fi) + Lsyriwd) - J 
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and Fi{p,e), Hi{p,e) are suitable continuous functions. 

With this notation in mind, we state the next proposition which claims that a large value 
of k succeeds in rendering bounded the trajectories of the switched nonlinear system (PHj) 
provided that the sampling interval T is sufficiently large. 

Proposition 2 Consider system \21\] with initial conditions inPxWxZxYxE. Suppose 
assumptions (A0)-(A3) hold. Let k, be chosen as indicated in LemmaU\ Then there exist 
T* > and k* > such that for all sampling intervals T > T* and all k > k* the trajectories 
are bounded in positive time. 



Proof. System (PHj) is apparently identical to the closed- loop system already studied in [2], 
but with the exception that, every T units of time, is being "reset" as 

Wd{kT) = WiikT-) + 

and e is being "reset" as 

e{kT) = e{kT~) - yr{w^{kT-) + w^ik)^) + y^iw^ikT-)) . 

With this in mind, we proceed to study the behavior of (|23|) on the time interval [0,T). 
Observe that system 

P = Foip) (24) 

coincides with system fjl7|) . the only difference being that the component w of z is now 
written as and ^ is now written as ^. Thus, from Lemma ^ it can be asserted that in 
this system all trajectories with initial conditions in Z x H are attracted by the compact 
invariant set 

^={(z,Og AxM'^:e = r(z)}. (25) 

Moreover, by construction, the function Ho{p) vanishes on this set. 

Let V denote the domain of attraction of A. Then, as shown for instance in |3], for 
any given arbitrarily small e, it is possible to claim the existence of a continuous function 
V : V ^ W having the following properties: 

(a) V{p) = a p ^ A and V{p) > everywhere else, 

(b) V{-) is proper on V, 

(c) for some positive b 

V{p)<b ^ \p\A<e, 
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(d) for some positive g < b, V{-) is locally Lipschitz on the set Vg = {p E V : V{p) > g}, 

(e) for any p E Vg , 

limsup hv{(f){h,p)) - V{p)] < -1 
in which (f){t,p) denotes the flow of (j21I). 

Let now e < 62 (with ^2 defined after Lemma Q), pick the function V accordingly and note 
that by construction 

V{p) <b \p\a <e<d2 peZxE. (26) 

Moreover pick numbers a and 61, with a > b such that Z x H C ^"^([0, a]) (which is possible 
as V is proper on V), and g < bi < b, and set 

S = {peV : bi< V{p) < a + 1} . 



Let c and A be positive numbers such that 

\y — Uriw)] < c for all y eY and w G Wq 

and 

IVriw) — yr{wd)\ < A for all Wd eW and w G Wq 

and note that 

\e{0)\ = \y{0) - y^iw^mi = \y{0) - y,{w{0)) + y,{w{0)) - y^iw^mi 

< \y{0) - yr{w{0))\ + \yr{w{0)) - 2/r(^d(0))| < c + A . 

Finally, let / and h be defined as 

/ := max |Fi(p, e)| h := max \Hq{p) + Hi{p, e)e\ . 

pev-i([o,a+il) pev-i([o,a+ii) 

|e|<c+A+l |e|<c+A+l 

Claim 1: There exists T* > and k* > such that for any T >T* , any k > k* and any 
£>0 

\e{iT)\<c + A m+l)T-)\<c ^ 

p{iT) eZxE p{{i + i)T-) e Z X E ■ ^ ^ 

To prove the claim, we derive two basic inequalities. The first one is about the dynamics 
of e. Looking at the bottom equation of (|23|) and bearing in mind the definition of h, standard 
arguments can be used to claim that - if on some time interval [to, (^ + ^)T), to ^ ^T, the 
state of ^ satisfies \e{t)\ < c + A + 1 and p{t) G V-^{[0, a+1])- then 

\e{t)\ < e-'^(*-*°)|e(to)| + ^ Vt G [to, {i + 1)T) . (28) 
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The second inequality is about the dynamics of p. Looking at the top equation of (j23|) 
and taking the Dini's derivative of V{p{t)) we see that - if on some time interval [to, (^+1)^), 
with to > iT, the state of satisfies \e(t)\ < c + A + 1 and p(t) E S - then 

D-'Vipit)) < -1 + Lvfmi < -1 + Lvfie-'^'-'°^\eito)\ + ^) Vt G [to, (i + l)T) , 

(29) 

in which Ly is such that \V{p) — V{q)\ < Lv\p — q\ W p,q E S. From this, using the 
appropriate comparison lemma, we obtain 

V{p{t)) < Vipito)) - (1 - ^^){t - to) + MMM vt e [to, {i + 1)T) . (30) 

In order to be able to prove Claim 1, we need this auxiliary result. 
Claim 2: there exists k** such that for any k > k**, T >T* and i>0 

\e(m\<c + A |e(t)|<c + A + l ^ , 

(31) 

This claim can be proved by contradiction. As a matter of fact suppose that ()3H1 is not true, 
namely that there exists a time T' e [iT, {i + 1)T) such that either |e(T')| > c + A + 1 or 
V{T') > a + 1 (with a mild abuse of notation we write V{-) for V{p{-))). By continuity of 
the trajectories with respect to time, this means that there exists a time T" > £T, T" < T' , 
such that either 

|e(t)|<c + A + l te[£T,T") |e(T")| = c + A + 1 and V{t)<a + l t e [iT,T"] 

(32) 

or 

V{t)<a+1 te[iT,T") V{T")=a + l and |e(t)|<c + A + l t e [iT,T"] . 

(33) 

But if (j32|l were true, by (j28|) taking to = and k > 2h, we would have that 



|e(T")| < e-'=(^"-^^)|e(£T)| + ^ < c + A + ^ < c + A + I 

K rh jL, 

which contradicts |e(T")| = c + A + 1. A similar contradiction would be obtained if (jHHj) 
were true. In fact, let to, < to < T" , be any time such that V{to) = a and V{t) > a for 
all t G [to,T"]. Using (jHUI) and taking k > max{Ly//i, 2Lv//(c+ 1)}, we have that 

VipiT")) < a - (1 - ^)(T" - to) + M|±l) < , + 1 

which contradicts V{T") = a + 1. From this. Claim 2 follows by taking 

k** = max{2h, Lvfh, 2Ly/(c + A + 1)} . 
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o(End proof Claim 2) 



Having proven that \e{t)\ < c + A + 1 and p{t) G l^~^([0, a + 1]) on the entire time interval 
[iT, {£ + 1)T), the desired Claim 1 follows again from arguments based on the bounds (j2HI) 
and ()30|) . To show that a large value of k succeeds in recovering \e{{£ + 1)T~)\ < c, just set 

,2h 1 2(c + A), 

k > maxj — , — In 1 

c T c 

in the estimate ()28p evaluated with = ^T. 

To show that a large value of k succeeds in recovering p{{£ + 1)T~) G Z x H we prove 
that if 

T >T* = 2{a-b) + l 

then V{p{{i + 1)T-)) < b which, by implies + 1)T~) G Z x S. To this end, we 
distinguish two cases. If p(t) G S for all t G [£T, (£ + 1)T), just set 

> max{2Lvf{c + A + 1) , 2Lvfh} 

in the estimate (jHUj) with to = to obtain V{{£ + 1)T~) < 6. The condition p{t) G 5 
for all t G [£T, (£ + 1)T) can be violated if and only if there are times in [iT, {£ + 1)T) 
at which V{t) < hi (recall the definition of S and the fact that V{t) < a + 1 for all 
t G [IT, {i + 1)T)). Let this be the case. If V{{£ + 1)T~) < h < b, the claim is triv- 
ially true. If not, let T' G [iT, {i + l)r) be such that V{T') = hi and bi < V{t) < a + 1 for 
all t G [T', {i + 1)T). On this time interval one can still use ()30|) with to = T' to conclude 
that iik> Lyfic + A + - bi) then V{t -T')<b for all t G [V, {i + 1)T) from which 
the Claim 1 follows. o(End proof of Claim 1) 

With this result at hand. Proposition 2 can be easily proved by subsequent applications of 
Claim 2. To this end note that p(0) G Z x H and, by definition of c and A, |e(0)| < c < c+ A, 
from which Claim 1 evaluated for £ = yields |e(T~)| < c and p{T~) G Z x S. Suppose now 
that, for some i > 0, \e{iT~)\ < c and p{iT~) G Z x H. Then, after the switch, 

\e{iT)\ = \y{iT) - yrMiT))\ = \y{iT~) - yrMiT))\ 

= \y{iT") - yr{wd{iT~)) + yr{w^{iT~)) - yr{w^{iT)\ 

< \e{iT-)\ + \y,{wa{iT~))-y,{wa{iT))\ < \e{iT-)\ + A < c + A 

and, by item (i) of Proposition^ Wd{iT) G W. The latter, by bearing in mind the definition 
of the set Z and the definition of p, yields p{iT) G Z x S from which the result of Proposition 
12] follows. < 

Proposition 2 shows that trajectories of the controlled system remain bounded if the 
time interval T exceeds a minimum number T* (minimal "dwell-time" ) which depends on 
the parameters of the controlled system and on the sets of initial conditions. This, however. 
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may be in contrast with relation (jlUj) which, bearing in mind that M(-) is an increasing 
function, requires that the samphng interval T is small enough. In the case the dwell-time 
T* is not compatible with (fTUI) . a simple modification of the decoder structure helps solving 
the problem. Let Nh be given, let T denote the minimal value of T compatible with (fTUI) 
and let i be any positive integer such that 

iT>T* . 

Consider now a "second level" decoder dynamics defined as 

< = s{w',) w'M = wM (34) 

whose state w'^ is periodically reset, every iT units of time, to the value of the "first level" 
decoder dynamics ((HI), that is as 

w'^ikiT) = Wd{k£f) for all A; > . 

In other words w'^ provides an under-sampled version of the first-level decoder dynamics (jH)) 
with the under-sampling period £T such that the constraint on the minimal dwell-time is 
respected. 

The same arguments used to prove Proposition 1 show that w'^(t) e W for alH > and 
furthermore 

lim \w'^{t) -w{t) \ = 

t^oo 

with uniform convergence rate. 

Consider now the regulator ()2Up in which e is defined as in p9|) but with Wd replaced by 
w'^. It is easy to realize that all the analysis carried out in this section can now be repeated 
by replacing the decoder dynamics (jHI) by (j34p . the decoder variable by w'^, and the time 
interval T by iT. In particular by Proposition 2 and by the fact that iT > T* we conclude 
that the trajectories of the controlled system are bounded if k is chosen sufficiently large. 
This and the fact that w'^(t) converges asymptotically to w(t) with uniform convergence rate 
are the crucial properties needed to conclude that the proposed regulator solves the problem 
in question as precisely stated and proved in the next subsection. 

5.2 The tracking error converges to zero 

To prove that the tracking error converges to zero, it is useful to observe that, if the coordinate 
y of (PT|1 is replaced by 

e = y- Vriw) 
the system in question can be also rewritten as 

w = s{w) 

z = f{z,e + y,{w),n) 

i = ^,{^) + G{-ke) + G{-ke) ^ ' 

e = q{z,e + y^{w), ji) — LsVriw) + VE_ — ke — ke 
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having set 



e . 



The same change of variables used to put (j22j) in the form (|23|) yields now a system of 
the form 

V = Fo{p) + Fi{p,e)e 

e = Hq{p) + Hi(p, e)e — ke — ke , 

in which p = col(/x, w, z,^) and Fq{p), Fi{p, e), Hq{p), Hi{p, e) are the same as in This 
system can be viewed as a "perturbed" version of system 

p = Fo{p) + Fi{p,e)e 

e = Hq{p) + Hi{p, e)e — ke 

whose asymptotic properties have been investigated in jl]. 

The following result is a minor enhancement of the main result of jl]. Let V{-) be the 
positive definite function introduced in the proof of Lemma ^ set V = {p : V{p) < a} with 
a chosen so that PxWxZxSgV, and set £" = {e : |e| < c}. Moreover, let A be the set 
defined by 

Lemma 2 Consider system ^37}) in which Fq{p), Fi{p,e), Ho{p), Hi{p,e) are defined as be- 
fore and initial conditions are taken in V x E. Suppose assumptions (A0)-(A3) hold. Let 
K he chosen as indicated in LemmaU\ Then there is k* such that, if k > k* , the following 
holds: 

(i) the positive orbit ofVxE under the flow of (|J7| ) is bounded and 

lim|p(t)U = 0, lime(t) = 0. 

4— >oo t— >oo 

(ii) for any e > 0, there exist numbers Si > and ^2 > such that, if |poU ^ C'l^d 
\go\ < Si, for any continuous function u{t) satisfying \u{t)\ < 82 for all t > 0, the solution 
p{t), e{t) of the perturbed system 

p = Fo{p) + Fi{p,e)e , . 

e = Hoip) + Hi{p,e)e-ke + u{t) ^^""^ 

with initial conditions p{0) = po and e(0) = Cq satisfies 

\p{t)\A<e, \e{t)\<e, Wt > . 

Proof. Item (i) has already been proven in The proof of item (ii) consists in a minor 
modification of the arguments used in |3] to prove item (i). Consider the Locally Lipschitz 
Lyapunov function U{p) defined in which - for some ai < 1, A > and L > - satisfies 



ai\p\A<Uip)<\p\ 



A ■ 
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and 

D^U{pit))<-\U{pit)) + Lf\eit)\ 

along the integral curve p{t), so long as p(t) remains sufficiently close to A (see (24) in 
Consider now, for ()38p. the candidate Lyapunov function 



which trivially satisfies 



W{p,e) = ^{U\p) + e') 



Y(\P\A + e')<Wip,e)<li\p\^^ + e^). 



Taking its Dini derivative along the trajectories of (joHj), we can obtain an estimate of the 
form (we omit the argument t for convenience) 

D+W{p,e) < -XU^ip) + LfU{p)\e\ + P\pU\e\ - (A; - ^)|ep + |M||e| 

so long as p(t) remains sufficiently close to A, in which f] and k are fixed positive numbers. 
Using the above estimates for U{p), it is easy to deduce that 

Clearly, there is a value k* > and a number a > such that, if k > k* 

and this yields, using the appropriate comparison lemma, 

W{p{t),e{t)) < -e-'''W{po,eo) + ^ max \u\t)\ , 

for all t > 0. From this, the result follows by standard arguments. < 
We are now ready to prove the main result of the paper. 

Proposition 3 Consider system \21\) with initial conditions inPxWxZx'ExY. Suppose 
assumptions (A0)-(A3) hold. Let n he chosen as indicated in Lemma^and k as indicated 
in Proposition^^ and in Lemma\^ Then 

lim e(t) = . 

t— >oo 
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Proof. As already mentioned, system (|2ip can be written in the form (|36|). Moreover, if 
the initial condition of (j2ip is taken in the set P x W x Z x E x Y , the corresponding initial 
condition of (jHUj) is in P x E. In view of the result of Lemma |21 item (ii), the result is proven 
if we are able to show that, given any pair of numbers 5i and 62, there is a time t such that 

\p{i)\A<Si, m\<6, (39) 

and 

\ke{t)\<52, yt>i. (40) 
To this end, recall that by definition 

e = e-e = -yr{wd) + yr{w) . 

Hence, since ?/i.(-) is continuous, in view of Proposition Q we have that limj^oo e(t) = 0. As 
a consequence, there is a time t*, dependent on 62, such that (|4(jp is fulfilled for all i > t* 
(note that the coefficient k - which is possibly a large number - is now fixed). Thus the only 
critical issue is to make sure that ()39|) holds for some t > t*. To this end, one can use the 
following argument, suggested in 0. 
Consider a system 

x = f{x) + u{t) (41) 

in which /(■) is locally Lipschitz and u(t) is a piecewise-continuous function. Let x{t, tQ.xo, u) 
denote the integral curve passing through xq at time t = to. Suppose u{t) satisfies 

lim u{t) = , 

and that, for a given xq and a given T > 0, there is a compact set X such that 

x{t,0,XQ,u) e X Vt>0, (42) 

x{tJT,x{iT,0,xo,u),0) e X yt>iT, V£gN. (43) 

Claim 3: for any t2 > and any 6 > 0, there is a i* such that, for all i > £*, 

\x{t2 + iT,0,Xo,u)-x{t2 + iT,£T,x{iT,0,Xo,u),0)\<6. < (44) 

To prove this claim, set x{t) = x(t,0,Xo,u) and x{t) = x{t,iT,x{iT,0,xo,u),0). For all 
t > iT, we have 

x{t) = x{£T)+ I f{x{s))ds+ I u{s)ds 
Jer J IT 

and ^ 

x{t) = x{iT) + [ f{x{s))ds 

J IT 
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Since /(■) is locally Lipschitz and X is compact, there is L such that \f{x) — f{ii)\ < L\x — y\ 
for all X and y m X. Then 

\x{t) - x{t)\ < L ! \x{s) - x{s)\ds + f \u{s)\ds. 
Jet JiT 

Since u(t) converges to as t — > oo, given any e there exists Tg such that \u{t)\ < e for all 
t > Ti;. For any iT > we have 

\x{t) - x{t) \<l[ \x{s) - x{s)\ds + e{t - iT) , 



from which Gronwall-Bellman's lemma yields 



with t > iT, hence 



\x{t)-x{t)\ < i(e^(*-^^) - 1) 

1j 



\x{iT + t2) - x{iT + t2)\< y (e^*2 - 1) 



Given any ta > and 5 > 0, set e = SL/{e^^^ - 1) and choose i* so that i*T > T^. This 
proves Claim 3. 

With this result at hand, it is easy to complete the proof of the Proposition. To this end, 
we set X = col(p, e), and identify system ()37|) with a system of the form 

X = fix) (45) 

and system ()36p with a system of the form (PT|) . The assumptions under which the previous 
claim holds are satisfied, with X taken as the set of all [p, e) such that p G V^~^([0, a + 1]) and 
|e| < c + 1. In fact, the proof of Proposition |21 shows that trajectories of (|36p are contained 
in this set X, i.e. that condition (|42p holds. The same proof also shows that at each time 
£T, where T is the sampling interval (of the over-sampling period according to the discussion 
at the end of the previous subsection), trajectories of (jHBj) are such that p G y~^([0, a]) and 
|e| < c. This and the fact that system coincides with (jHUj) with n = in turn guarantees, 
by Proposition 121 that also condition ()43p holds. Since x{t,0,x{iT,0,XQ,u),0) is a solution 
of (j37|) with initial condition inVxE, we know from Lemma |2l item (i), that - given any 
6i > - there exists t2 such that the p and e components of x{t2, 0, x{iT, 0, Xq, u), 0) satisfy 

|P(^2)U<|, |e(t2)|<|. 
Using now the Claim 3 with 6 = 6i/2 and the fact that 

x{t2 + iT, iT, x{iT, 0, xo, u),0)= x{t2, 0, x{iT, 0, xq, u),0) , 
we deduce that the p and e components of x{t2 + iT, 0,xo,u) satisfy, for all i > i* 

\p{t2 + iT)\A<Su \e{t2 + iT)\<6i. 
This is what was needed to complete the proof of the Proposition. < 
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6 Simulation Results 



We consider the problem of synchronizing two oscillators located at remote places through 
a constrained communication channel. The master oscillator (playing the role of exosystem) 
is a Van der Pol oscillator described by 

wi = W2 + e{wi + awf) 

W2 = —Wi 

whose output Href = W2 must be replied by the output ?/ of a remote system of the form 

y = u. (47) 

Simple computations show that, in this specific case, the steady state control input u^s 
coincides with u^^ = —wi and the assumption (A3) is satisfied by 

ii = 6 
6 = /(6,e2) 

Uss = 6 

where /(^i,6) = "6 - ^{^2 - •ia^K^) through the map 

t{w) = ( —Wi —W2 + e {wi — awf) )^ 

We consider a Van der Pol oscillator with e = 1.5 and a = 1. The regulator (j2(jp is tuned 
choosing k = 3, G = (12 36) and k = 8. We consider two different simulative scenarios 
which differ for the severity of the communication channel constraint. In the first case we 
suppose that the number of available bits is Nb = 2 yielding, according to © and to the fact 
that r = 2, = 2. In this case, for a certain set of initial conditions, condition © is fulfilled 
with T = 0.15 s. In the second case the available number of bits is assumed A^f, = 4 from 
which © and Q yield a bigger A^ and T respectively equal to A^ = 4 and T = 0.5 s. The 
simulation results, obtained assuming the exosystem ()46|) and the system ()47|) respectively 
at the initial conditions w{0) = (1,0) and y{0) = 5, are shown in the figures for the 
first scenario and figures 12113 for the second one. In particular figure Q (respectively shows 
the quantized variable Wg transmitted from the encoder to the decoder and used to reset 
the respective dynamics according to the rule described in Section 01 Note that in the first 
scenario represented in figure ^ each of the two components of the vector Wq, taking value 
in the set { — 1/2, 1/2}, can be transmitted using 1 bit. On the other hand in the second 
scenario, represented in figure El the transmission of each component of Wq, whose value 
are in the set {—3/2, —1/2, 1/2, 3/2} requires 2 bits. Figure El (respectively E]) shows in 
the left-half side the error, as a function of time, between the exosystem and the encoder 
(decoder) state and in the right-half side the phase portrait of the Van der Pol oscillator with 
overlapped the actual state trajectory of the encoder (decoder). Finally figure El (respectively 
inj plots the the tracking error e(t) = y{t) — y^-it) on the left-half side and the control input 
u(t) on the right-half side from which it is possible to see that in both the control scenarios 
the synchronization between (|46|) and (|47p is achieved. 
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Figure 1: First control scenario (A^ = 2, T = 0.15 s): behavior of the encoded variables 
Wgi{t) (left) and Wq2{f) (right). 



7 Conclusions 

We have discussed the problem of asymptotically tracking a reference signal which is gen- 
erated by a remote exosystem, and transmitted through a finite bandwidth communication 
channel. Although only an estimate of the actual tracking error is available to the regulator, 
a suitable choice of the controller parameters allows us to achieve the control goal while 
fulfilling the constraint on the bandwidth of the channel. 
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Figure 2: Second control scenario (A^ = 4, T = 0.5 s): behavior of the encoded variables 
Wgi{t) (left) and Wg2(t) (right). 
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Figure 3: First control scenario (A^ = 2, T = 0.15 s). Left: time behavior of w{t) — Wc(t) 
{w(t) —Wd(t)). Right: phase portrait of the exosystem (dotted hne) and trajectory (^01,^02) 
(sohd hne). 
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Figure 4: Second control scenario (A^ = 4, T = 0.5 s). Left: time behavior of w(t) — wdt) 
{w(t) —Wd{t)). Right: phase portrait of the exosystem (dotted line) and trajectory {wei,We2) 
(solid hne). 
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Figure 5: First control scenario (A^ = 2, T = 0.15 s). Left: time behavior of the reference 
trajectory yr{t) (dotted hne) and of the controlled output y(t) (solid line). Right: time 
behavior of the control input u(t). 
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Figure 6: Second control scenario (A^ = 4, T = 0.5 s). Left: time behavior of the reference 
trajectory |/r(t) (dotted line) and of the controlled output y{t) (solid line). Right: time 
behavior of the control input u(t). 
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